Up to and beyond ninth order in opacity: Radiative energy loss with GLV 
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A new examination of the GLV all-orders opacity result for radiative energy loss is presented. The 
opacity expansion is shown to be a Dyson expansion of a Schrodinger-like (or diffusion) equation, 
a form also found in BDMPS-Z-ASW, AMY and Higher Twist formalisms of radiative energy loss. 
A new, efficient numerical evaluation of the GLV all-orders result gives results up to and beyond 
ninth order in opacity. This presents a solution to this Schrodinger-like equation to almost arbitrary 
accuracy for lengths and densities applicable at RHIC and LHC. The first order is found to be 
most important, although higher orders implemented in a geometry integration will be needed for 
quantitative jet tomography. 

PACS numbers: 12.38.Mh; 24.85.-(-p; 25.75.-q 



I. INTRODUCTION 

Jet quenching theory was buoyed by early success at 
RHIC. This was based on the application of calcula- 
tions based on perturbative quantum chromodynamics 
(pQCD) to high momentum particles (jets) traversing a 
quark-gluon plasma (QGP) at mid-rapidity produced in 
Au-Au collisions at ^i = 62 - 200AGeV % UMM- 
Reasonable choices of input parameters could predict or 
explain high momentum pion and eta Raa{pt) data out 
topT«20GeV 0,0,0]. 

These calculations were based on radiative energy loss. 
Coherence effects due to the soft emitted gluon make ra- 
diative calculations a difficult business. There are cur- 
rently four different calculations on the market for ra- 
diative energy loss in a QGP: (^)GL\Mi, Jl [M [ll 111, 
BDMPS(-Z-ASW) P^[M[ll4ilJp,[il[l^^ 
m. Higher Twist [2l,l2El2ll23rAMY[2in2a & 
In spite of the similar roots of all the calculations, lit- 
tle agreement between groups can be found, in results 
or even use of language. The long term disagreement be- 
tween four similar schemes presents significant difficulties 
for the interpretation of experimental data at RHIC and 
the upcoming LHC. In the PHENIX study [sj, this is 
represented as an approximate factor ten uncertainty in 
the extracted medium parameter. Clearly we have not 
yet achieved quantitative jet tomography. 

The strongly coupled techniques of AdS/CFT have 
rece ntly surged onto the jet quenching scene [33, Isj, 
KMI . ISa . ISq . Many of these implicitly reject the no- 
tion of radiative energy loss, the emission of relatively 
soft quasi-particles. The heavy quark puzzle, the un- 
expectedly [37| strong quenching of heavy quarks seen 
throu gh t he low Raa{vt) of single non-photonic elec- 
trons [sl, [si, SO, El SI |43|, |4||, has also motivated 
broader investigation of energy loss mechanisms. Pos- 
sible 'solutions' include a re-examination of coUisional 
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energy loss ka |46|, |47|, |48|, heavy meson dissociation 
mechanisms [4!|| , resonances in the medium J 50l . ISlL l52l | , 
or the above-mentioned AdS/CFT techniques. A con- 
sistent theory of radiative energy loss is important to 
indicate whether further mechanisms are necessary. 

This paper begins to bridge the gap between radiative 
energy loss formalisms. This begins in section |lT] by em- 
phasizing the common foundations of all the calculations 
of gluon radiation. The GLV result is expressed as an 
opacity expansion. This is shown to be a Dyson expan- 
sion of a Schrodinger-like (or diffusion) equation. This 
Schrodinger-like equation is not confined to thick plas- 
mas or uncorrelated medium assumptions, and has been 
used in other radiative energy loss formalisms. 

This in itself is not new. But the use of the opacity ex- 
pansion is twofold: in the derivation of the Schrodinger- 
like equation (including the explicit demonstration of its 
color triviality), and as a means to its numerical solution. 
The GLV result was expressed in closed form to all orders 
in opacity. A new numerical evaluation of this result is 
presented in section IIIIl which is able to find results up 
to and beyond ninth order in opacity. 

The numerical results presented here use a Gyulassy- 
Wang (GW) static color center model of the medium. 
Possible improvements to this model have been sug- 
gested Djordjevic and Heinz's dynamical medium calcu- 
lations [53|, |5J]. The formalism and numerical solution 
presented here are not limited to this GW model, and 
future studies can explore further possibilities. For ex- 
ample, Vitev has given all-orders expressions for initial 
state, final state and infinite time radiation |55l |. 

In this way, for lengths and densities applicable at 
RHIC and LHC, all the significant terms in the series 
are found. The interference with the creation radiation 
is crucial to an accurate evaluation of the radiation spec- 
trum. Other numerical solutions to similar Schrodinger- 
like equations have been presented by other groups, and 
in future these solutions can be compared to the method 
of solution found here. There is no need to sit still with 
a fixed thin or thick plasma approximation. Once this is 
solved for a fixed length and density, the incorporation 



of these results into geometry integrations can mark the 
beginning of the era of quantitative jet tomography. 



II. THE GLV OPACITY EXPANSION 

In GLV II [i3| a reaction operator was derived, used 
to iterate from order to order in opacity and to prove 
the 'color triviality' of this procedure, and to derive a 
closed form expression for an all orders in opacity result. 
Here, we look again at the operator recursion in opacity, 
show its equivalence to a Schrodinger-like equation (for 
both thin and thick plasmas, correlated and uncorrelated 
media) similar to those used in other radiative energy loss 
formalisms, and discuss the significance of this result. 



A. The Schrodinger-like equation 

The GLV recursion relation is derived in GLV II [lOJ . 
This section continues on from that derivation. It uses 
eikonal kinematics for both the jet and the emitted gluon 



k±-^uj , E:^ q,q° 



(1) 



and a small x approximation - the energy of the emitted 
gluon is much less than the energy of the jet. The lon- 
gitudinal momentum of the gluon at different stages is 
written as oj,-. 



Wo 



2xE 



, uJi 
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2xE 



(2) 



This assumes that the energy transfer with the medium 
is zero, or at least negligible compared to the other terms. 
The GLV result expresses the relation between succes- 
sive orders in opacity: 



P„(k) = Ca [Pn-i{k - q„) - P„-i(k)] 
- 2CaB„. (Ree-'""^"e^'i-^I„_i 

+ Sn,lCACR\Bif (3) 

where Pri(k) is, up to multiplicative constants, the emis- 
sion probability of a gluon of momentum k (where k is 
momentum transverse to the direction of the jet). 

The first line is identified as a classical cascade. Once 
the emitted gluon has decohered from the jet, it collides 
with the medium and diffuses in 2 dimensional k {—k±) 
space. This does not change the total number, it inte- 
grates to zero. 

The second line is identified as a quantum cascade. It 
acts as a source term into the classical cascade. Each 
possible collision evolves the amplitude and its complex 
conjugate, giving evolution in I„. This recursion is 

I„ = Ca (e»(-o--„)^„e*'i-^ - l) I„_i 

- J„_iCACflBie^"«^i (4) 



The first line iterates from collision to collision, the sec- 
ond only applies to the first collision. The source from 
the second collision is smaller than (LPM suppression) 
or equal to (the incoherent limit) the first collision. 

The last line in Eq. ([3]) is a source from the first col- 
lision. We rewrite Eq. (I3|) as a classical cascade with a 
source term from a quantum process, 

Pn(k) = -P„xlassical(^) 

+ '-'n, quantum''*) + "n,i^i, quantum C*) (5) 

'5'n quantum ^^^^ ^'^ written as 5„. The iterations need 
initial conditions, ie Pq and Iq. These differ depending 
on whether the incoming jet is on-shell, ie created at 
t = — oo, or created at t — to. For the on-shell jet, 

Po(k) = , lo = (6) 

For the finite time case, 

Po(k) = CflH2 , lo = -CflHe'"«^« (7) 

Note how a Ca is always associated with an interaction 
with the medium, the Cr is associated with the emission 
vertex. 

B. Rearranging the formulae 

The shift operator appears in the iteration for I„, 
Eq. (U). This form will appear frequently: 



O^^CAie 



i{uja- 



"'e 



iq^.b 



1) 



(8) 



The first term corresponds to the phase shift from the 
(soft) gluon colliding with the medium. The second term 
corresponds to no phase shift from the (hard) jet col- 
liding with the medium. At this stage, the operator Oi 
is a function of the position z^, the momentum trans- 
fer Qi and the gluon momentum k. Averaging over q^ 
and Zi will come later. The shift operator is remarkably 
simple, hiding its complex derivation through the evalu- 
ation of direct and virtual QCD diagrams. This form is 
specific for the eikonal approximation, with no (or negli- 
gible) energy exchange with the medium. Note that the 
shift operator at Zi is only dependent on the momentum 
transfer q^, and independent of the rest. 

First, we want to simplify the beginning and end of 
the series giving Sn- B; can be rewritten with use of the 
shift operator on a hard vertex: 



CaB^c" 



-0,He" 



(9) 



In this way, we can express the first and last collisions in 
terms of the shift operator 



J/0„ 



Ii 



"e 



icin.h 



(10) 

(11) 



where 



gives 



J/(z„) = -CrUc-^'^"'" (12) 

J,(zo,zi) = -CKHe-^-"-i 

for a jet created at t=— oo (13) 

J,(zo, z,) = -CflH(e-*"«^i - e-^'^"^") 

for a jet created at t=0 (14) 

Now we have the quantum cascade contribution to the 
gluon spectrum in a simpler form: 

Sn{k) = -2J/(Z„) OnOn-l . • . Ol J.(zo, Zi) (15) 

This wiU be averaged over g^ and Zi. First, we look at 
Qi, and define 

A(zi) 



A(Zi) 



(16) 



This is the momentum averaged shift operator. It is lo- 
cal in z. The potential has been normalized to give V, 
with the mean free path beneath it (the Ca has been 
absorbed into this A(z)). Remembering that (...) rep- 
resents averaging over momentum transfers (not the z^s 
yet), the expression for Sn now becomes 

(5„(k)) --2Jy-(z„)(0„)(0„-i)...(6i)J.(zo,zi) (17) 
We now integrate over the ZiS. 

('5„(k))q^^^; ^-'2 dzi dZ2... dznX 

J Zq J Zi J Zn 

J/(z„)(d„)(0„-i)...(Oi)J»(zo,zi) (18) 

All the Zi integrations pair up with their Oi, except the 
first and last, zi and z„ - as the currents Ji(zo,zi) and 
3 f{zn) introduce additional z dependence. Now the aim 
is to pull the dzn integration out to the front. Using 

dz„^i / dzn = / dzn dzn^i (19) 

-Sti — 2 '^ ^n — 2 '^ ^n — 2 '^Zn — 2 

this gives 

{S.a{k))q.^,_. = -2 / dz„J/(z„)(0„)x 

J Zo 



dzi / dz2 ■ ■ ■ I dZn-lX 

J Zi J Zn-2 

(O„_i)...(6i)j,(zo,zi) (20) 



Defining ij}^ 



i)^Xzo,Zn)^ / dzi dz2... dZn-lX 

J Zo J zi J Zn-2 

(O„_i)...(6i)j,(zo,zi) (21) 



(5„(k))g,,,, =-2 / dz„J/(z„)(O„)'0„(zo,Zn) (22) 



C. Summing over opacities 

The result will be summed over all possible opacities, 
from n = 1 ^ oo. The relevant opacity weights come 
out from the calculation, as the A^s have already been 
included. The result is 

CO CO -.jg^ 

^(S'„(k))5,_2, = -2^ / dz„J/(z„)(6„)t/>„(zo,z„) 

n=l Ti=l'^2o 

(23) 
The Zn variable as a dummy variable can be replaced by 
Ze, and pulled out of the sum, giving 

^(5„(k)),„,, = -2 / dZeJ/(Ze)(Oe) 5]tA„(zo,Ze) 
n=l -^^0 n=l 

(24) 
Defining the full wave function 

oc 

■0(zo,Ze) = X!^«(^"'^«) (^^) 

n=l 

gives the equation for ip as 

/•-Se fZ^ 

dZn-lX 

n-2 

(O„_i)...(0i)j,(zo,zi) (26) 

Written out as a series, this becomes 

pZe 

■j/'(zo,Ze) = / dzi . .. (Ol)Ji(zo,Zl) 

J Zo 

+ dzi dz2(02)(Oi)Ji(zo,zi) 






+ / dzi dz2 dz3{03){02){Oi)3,{za,zi) 

J Zo J Zi J Z2 

+ ... (27) 

This is familiar from scattering theory, reminiscent of a 
Dyson expansion. 



Xp{zo,Ze) =T 



exp 



dz{d){z) 



(28) 



Hence, ip{zQ, Zg) satisfies a Schrodinger-like equation 

|-tA(z,k)-(0)(z,k)tA(z,k) (29) 

oz 



D. The Hamiltonian in impact parameter space 

In q-space, the form of the operator is unfamihar, but 
a transformation into impact parameter space makes the 
Schrodinger analogy clearer. 

A change of phase to 



rpiz,k) = e''^°H(f>{z,k) 



gives 



(30) 

(31) 

(32) 
(33) 



Vfc+^T^(^(^'b)-l)U, (34) 



e'^^-^ipiz, k) = ie'"'^</)(z, k - q) 
and a Fourier transform 

^(z,b) = f (fke-'^-^4>{z,k) 

V{z,h) = Jd\e-^'^-^V\z,q) 

gives the equation 
.90 f -1 



' dz \2xE ''' ' "A(z) 



This is a Schrodinger-like or diffusion equation in two 
dimensions with an imaginary (ie absorbing) potential. 
The Hamiltonian is 



H = 



-1 



1 



2xE 



A(z) 



iV{z,h)-l) 



(35) 



This was expressed by Zakharov [56| as 



rjf \ 1 9'^ . n(z)g3(p,z) 

^^") = -2M(^a^-^ 2 ('') 

where p = b, and a^ is the 'cross section for interaction 
of the qqg system with a scattering center'. This can be 
used in path integral form, with the Lagrangian 

C=^M{x)p^~^zn{z)a3ip,z) (37) 

Wiedemann describes this eikonal expression as the non- 
Abelian Furry approximation [l9|. The p^ kinetic term 
comes from the eikonal u/q = k'^/2xE - terms beyond 
the eikonal will spoil the simple diffusion equation. The 
AMY result [20| is expressed as an integral equation, 
which can be recast as a Schrodinger-like differential 
equation. The Higher Twist expansion has been 're- 
summed' recentl y by Majumder, in application to photon 
Bremsstrahlung [53| , to give a similar diffusion equation. 



E. Interpreting the Schrodinger equation 

For the intuitive interpretation of the Schrodinger 
equation, it must be remembered that it is not simply 
a Schrodinger equation for the propagation of the jet. 
Zakharov goes into much detail of interpretation of this 



Schrodinger-like equation. From the GLV recursion, it is 
clear that the evolution is not just in the amplitude, but 
also in its conjugate. This is far simpler due to the color 
triviality produced - as shown in GLV, the reaction op- 
erator is color trivial, and hence the evolution considered 
simultaneously in the amplitude and its conjugate does 
not have complicated color dependence. 

An opacity expansion is a perturbative-like Dyson 
expansion of the full result. Hence, the GLV opac- 
ity expansion recursion is equivalent in form to other 
approaches that started from the Schrodinger-like ap- 
proach, ie BDMPS and Zakharov. The Schrodinger form 
does not assume a thick plasma, only the original approx- 
imated solutions of the equation do so - in a similar way 
to how a Schrodinger equation in non-relativistic quan- 
tum mechanics can describe one, two, or an arbitrary 
number of interactions. BDMPS derives the form of the 
equation, starting with a large Nc result and then giv- 
ing the general result Il3|. Zakharov uses path integrals 
to get a similar result J5y|. Similarly, AMY derives an 
integral equation form of their result [2^. Majumder de- 
rives a similar form for photon Bremsstrahlung in the 
higher twist expansion f57| . All of these approaches are 
equivalent at their basic level, but they differ in: 

- The arbitrary form of V(z,q), 

- The approximations made to derive the explicit 
form of the operator 0{z), 

- The method of solution of the Schrodinger-like 
equation, 

- The initial condition (finite time or infinite time) 
to iterate / evolve, 

- The implementation of fluctuations in gluon num- 
ber. 

On top of this, there are all the variations in numerical 
implementation and the treatment of the geometry of 
nuclear collisions. 



F. Solving the Schrodinger equation 

This section has reproduced and repackaged old re- 
sults, expanding upon (for example) one line of BDMPS 
which notes a Bethe-Salpeter form of an equation \13^ . 
The use of the derivation is in the solution to the 
Schrodinger equation. The GLV II result gives an all or- 
ders in opacity result in closed form. To get the result to 
arbitrary order, there is no need to go back to the Green's 
function or to diagrammatic methods, as done for exam- 
ple by Wiedemann [20] . The result is expressed in a form 
that can be implemented on computer, in matrices that 
can be expanded in computer memory as opposed to on 
paper. 

At first, in 2000 when the GLV II paper was written, 
the all orders result was of limited practical use. First 



order was quick, second order difficult, third order almost 
out of reach. Now, with the advances in computing and 
some numerical optimization of the evaluation, it is fea- 
sible to evaluate the result up to and beyond ninth order. 
Hence, solution to the Schrodinger equation is possible to 
almost arbitrary precision. Once a form of the equation 
is settled upon, this method of solution can be compared 
to other methods, such as an explicit evolution of the 
quantum wave function with discrete time steps. In this 
way, the differences between the current set of radiative 
energy loss formalisms can be resolved. 

The numerical solution to the opacity expansion is 
given in section IIIII The qualitative conclusion is sim- 
ilar to that for GLV II: for lengths applicable at RHIC 
and LHC, the first order is the largest contribution, the 
higher orders are smaller and cancel each other. This is 
perhaps surprising: similar to a perturbative expansion, 
we only expect quick convergence if the driving coupling 
is small. In this case L/\ is not small, it is in the range 
1^6. However, there is an additional way in which 
the second order can be much smaller than the first, one 
that is not available in introductory quantum mechanics 
courses. Our potential is imaginary, absorbing. If it is 
strongly absorbing, then the second order is much smaller 
than the first, even though L/\ is large. 

Zakharov noted a similar effect, when radiation was 
in the 'diffusion' zone [50|. In this region, the BDMPS 
method of solution, which makes what seems to be a rea- 
sonable assumption that if L/X ^ 1 then we need to 
consider n ^ 1, in fact misses out the largest contribu- 
tion, the n = 1 in opacity. He then goes on to make a 
different numerical solution to the Schrodinger equation 
in a subsequent paper }58| . In future, different numerical 
solutions to the same equation can be compared. 



an approximation to the full 3D evolution. The eikonal 
property is applied to both the jet and the emitted gluon 
{k±^ <C Lo). Going beyond this approximation would be 
valuable. 

GLV uses a leading log evaluation of the initial condi- 
tion, ie the zeroth order in opacity emission. However, 
there is a significant contribution to the energy loss from 
the region in which the approximation k ^ xE does not 
apply. Relaxing this assumption will involve relaxing the 
eikonal, coUinear kinematics - remembering that these 
kinematics have been applied not just to the jet but also 
to the emitted gluon. AMY evaluates the energy loss for 
an asymptotic, incoming on-shell jet - to be evaluating 
the same equation, the evolution from a hard process as 
in GLV will be needed. 



III. EVALUATING THE GLV ALL ORDERS 
RESULT 

A. Summary of formalism 

We turn to the numerical evaluation of the GLV all 
orders in opacity result. The result in GLV II [l3| is ex- 
pressed for the case of a massless jet created at a finite 
time emitting massless gluons. The assumptions and ap- 
proximations in the derivation include eikonal kinemat- 
ics, small X energy loss. Djordjevic applied the same 
formalism but kept explicit mass terms for the jet and 
emitted gluon, giving the 'DGLV result. This is very 
similar to the 'GLV result, but with an effective shifted 
mass term: 



' -mlil 



Pdglv 
Pglv 



Arx- 



2^2 



2 -0 



(38) 
(39) 



G. The future: beyond eikonal, small x, static 
media 



GLV is evaluated in the extreme small x approxima- 
tion; BDMPS and Zakharov have tried to push past this 
region. This involves a quantum cascade of both the 
emitted gluon and the emitting jet in impact parameter 
space. Relaxing the x -^ 1 assumption in the same way 
for both formalisms will lead to the same Schrodinger 
equation for both (as long as the color triviality of the 
result remains). 

GLV is evaluated with a Gyulassy-Wang model of 
static color centers. Moving beyond this involves im- 
proved modeling of the V{q'^,q) medium potential, simi- 
lar to recent work by Djordjevic [53l . |54| or in the AMY 
formalism (for example, [2^). The relaxation of the 
g° = assumption may alter the form of the shift oper- 
ator; however, its remarkably simple form suggests that 
a simple form may also emerge from a more advanced 
calculation. In scattering in normal quantum mechanics, 
the eikonal form emerges out of a more general result: the 
two-dimensional evolution in impact parameter space is 



The final result is expressed in terms of the interference 
between formation time factors. The inverse of the for- 
mation time is the (longitudinal) momentum, 



^{m,...,n) 



/3^ + (k-q,„-...q„)^ 
2xE 



(40) 



n will always represent the final scatter, m will vary from 
the the first up to the final scatter. The notation is simi- 
lar to that in GLV II, with bold letters (q, k) representing 
2D vectors in the transverse plane. Note that the w„ used 
before is just a special case of this version, with m = n. 
The denominator is the positive light cone momentum of 
the emitted gluon (assuming k±_ <^ uj). 

The radiated spectrum is expressed in terms of 'Cas- 
cade' terms: 



^(ili2---im) 



(k-q^i -q^a - ...-qj„J 



/32 -h (k - q,i - qj2 - 



<l^J' 



(41) 



with l3 being either Pglv or Pj^glv from above. This 
represents the shifting of the momentum of the radiated 
gluon due to momentum kicks from the medium. We 



define the 'Hard' term as a special case of C without any 
momentum shifts: 



H 



fP + k2 



(42) 



The 'Gunion-Bertsch' scattering terms are linear combi- 
nations of C: 

with the special case 

B, = H - C, (44) 

The full result, differential with respect to x and k 
(which is uniform over ipk), for the radiation spectrum at 
nth order in opacity is expressed as 



d^(") C^a, 1 



n(^K(qo-^^(q. 



dxd'^'k TT^ n! / -*■-*- V ^o(0 

i—l ^ y ^ / 
~'^^il,--- ,n) ' 2_^ B(m+1,--- ,n)(m,--- ,n)>< 



m— 1 



cos^u;(fc_..._„)Azfc-cos^W(fc_...^„)Azfc 



fe=2 



k=l 



(45) 



where ^2 = is understood and |5i(qi)p is defined as 
the normalized distribution of momentum transfers from 
the i*^ scattering center, \g{i) is the mean free path of 
the emitted gluon, Cr is the color Casimir of the jet. 

Noting the discussion in section III Al about classical 
vs quantum cascades, the expression for the source term 
(the only term that changes the energy distribution) at 
all orders in opacity is: 



dx cPh TT^ 71 !/-'■-'- V Xa(i) 

i—l ^ y ^ / 

X ( — 2C(i^...^„) • B„x 



ln(^WM->'M] 



COS ^ uj(k^... ,„) Azfc - cos ^ w(fe^... ,„) Azfc 



fc=2 



fe=l 



(46) 



1. Smooth background density 

For evaluation in a smooth background density, we 
change 



s/n(<i'<..^w«,.)-^^«..)i) 



dzi ' 



dZn Yli ^^^' 



z„-i J j^i 



2 w^(q»)-<52(q,) 



\[z) 



(47) 



where A(z) is is the distance (and time) dependent mean 
free path. 



2. Static color center model 

The GLV II result is not confined to a specific model of 
the jet-medium interaction. The ^(g°) in the Gyulassy- 
Wang potential has been used in the derivation, but work 
by Djordjevic and Heinz suggests that this strong as- 
sumption may not be needed to give the same form of re- 
sult [54] (with a slightly modified potential). In an HTL 
calculation, the unscreened magnetic interaction sends 
the mean free path to zero, violating the assumption 
A ^ l//i. However, these soft interactions do not con- 
tribute to much of anything, except collision number. An 
ad-hoc magnetic screening fj,m,ag ~ fJ-D pushes the mean 
free path up to « ifm again, without strongly affecting 
the (Ap) or (ql). 

For this evaluation, we stay with the simpler Gyulassy- 
Wang static color center model. For a specific uj, k, 
the mean free path does not affect the formation time. 
Hence, if the opacity series peaks at n = 1 for the static 
model (from Zakharov, if Tf ^ L), this should persist to 
other models with shorter mean free paths. 

For the Yukawa screened interactions, the differential 
gluon cross section in the local density approximation is 



<7{zi,qi) = aei{zi)\vi{qi)\ 



2 d'^aeiizi) 



d'^qi 






(48) 



where /i(zi) is the local Debye mass (that may vary if the 
system expands). 

The total cross-section, 0-^1 is related to the (gluon) 
mean free path, A = l/{paei) where p is the local density, 
and is given by 



2At2 



3. Relating X to q 



(49) 



The average q kick per collision, for q,nax ~ v6TE ^ 
/i, given by this potential is (q) = 0, but with 



E 
(q^)per collision = M^ log j^ 

hence the average q per unit distance is 
p\ E 



(50) 



E' 



(q )per fm = ^ l^S TT ^ ^l^g — (51) 



A ° /i^ fjp 



We could present the result in terms of q (which here 

7^ {q']_)) by the substitution 



A p^ 



(52) 



and expressing in terms of this q, but we still have a 
dependence on fi. In general, q does not uniquely give 



the energy loss; you either have to specify both q and ^, 
or both A and ^. In a kinetic theory, A and ^ are related 
to the temperature. Here, we will set A = Ifni, ji = 0.5 
GeV. 



where 



N N 



(56) 



B. The uncorrelated medium 

In general, A(z) can vary with z. But even if we have 
a box of constant temperature, A(z) = A, the spacing of 
collisions is mutually dependent. Taking the example of 
second order in opacity calculation confined to a length 
L, if the first collision happens at 0.9L, the spacing to the 
second is limited to less than O.IL. 

A significant numerical simplification is achieved by re- 
moving this mutual dependence, and evaluating an 'un- 
correlated medium'. This can be expressed as |10 | 



P{zi, 



Li Le{n) 



n 



J^e-A2j/Le(«) 



(53) 



The spacings of collisions are independent. This can only 
be considered an approximation to a box profile, which 
is then only an approximation to the varying density 
profile. The full numerical result is not limited to this 
approximation, and future work should examine the dif- 
ference between the box approximation and the uncorre- 
lated medium approximation. 

Eq. ([55)) converts the oscillating formation physics fac- 
tors in the full result into simple Lorentzian factors 



/m m 

(ipcos^W(fc,...,„)Azfc = Re ]^ - 

k=j k=j 



1 



K—J 

(54) 
In order to fix Le{n), we require that {zk ~ zq) = kL/{n + 
1) for both geometries. This constrains Le{n) = L/{n + 



Monte Carlo evaluation 



In this way, an estimate of the result and a rough esti- 
mate of the uncertainty on the result is found. The use 
of quasi-random numbers can often improve the rate of 
convergence of Monte-Carlo integration, but this is not 
pursued further here. 

The key to the multi-order implementation is the nu- 
merical evaluation of vector dot products between C and 
B. Considerable effort is saved by iterating these numer- 
ically as opposed to by hand on paper. The results pre- 
sented here use only the most basic non-adaptive Monte- 
Carlo integrators. This study is presented as a starting 
point, to show that high orders in opacity are well within 
the reach of current computing power. 

Only simple fixed length results will be shown. The 
gluon mean free path is fixed at Xg = Ifm, the Debye 
mass at iio = 0.5GeV. The q± momentum integral is cut 
at V&TE, with T=0.25GeV. 



IV. ALL ORDERS RESULTS 

Results were presented in GLV II |10| for summation 
up to 1st, 2nd and 3rd orders in opacity. In appendix 
1X1 explicit comparisons to some of the GLV II results 
are shown. Good consistency is found for first and sec- 
ond orders, while the third order result for [l3| is a little 
smaller than that found here. Error bars on all results 
are shown, as per the uncertainty estimate from Eq. (j56p . 

Results plotted in this chapter are for the one-gluon 
distributions. In order to produce an RAAiPr) result, 
one gluon distributions are usually input into a Pois- 
son convolution before being combined with initial spec- 
tra. Many of the results are compared to the Bertsch- 
Gunion spectrum, the incoherent radiation spectrum 
with a Gyulassy-Wang model and a small 0.1 GeV screen- 
ing of the k_L = qx pole. 



The GLV result is ideally suited to Monte-Carlo evalu- 
ation. For high orders in opacity, it is a high dimensional 
integral (either 2n or 3n depending on whether the uncor- 
related medium approximation is used or not) and should 
not be dramatically peaked in one small region (although 
it is highly oscillatory). 

Monte Carlo integration proceeds by repeated evalua- 
tion of the integrand. By throwing dice as to the value 
of the variables (over which the integration is being eval- 
uated) and averaging the results, an approximation is 
made to the result of the integral. For large numbers of 
evaluations, the result should converge to the true value 
of the integral as l/^/n. The result is expressed as [59| 



fdV « V{f) ± V 



ip) - if)' 

N 



(55) 



A. Results - dN/dxdk 

1. p = 100 GeV 

The interpretation of the results is complicated by the 
effect of kinematic limits, so we start with results for 
light quarks with p7'=100GeV. The individual orders in 
opacity (up to n=9) are presented in appendix section [Bl 
for L=6fm, E^lOOGeV, for x = 0.025,0.05,0.2,0.5,0.8, 
for both the full result (including the classical cascading) 
and for only the quantum source term (the part that 
is important for dN/dx). Here, we concentrate on the 
answers summed up to a specific order. Unless otherwise 
stated, a result shown uses the full result including the 
classical cascade, ie Eq. P5)) . 



To summarize: the first order result gives a good ap- 
proximation to the aU orders result for large k±, short 
lengths, large x. 

For all the dN/dxdk results (shown in Fig. [9l Fig. \lU\ 
Fig. [Til Fig- EH) 7 the all orders (summed up to 9th order) 
result deviates from the first order for small k±_ , and then 
follows the first order result for large k±. The point of 
reattachment depends on x, the energy of the emitted 
gluon. This can be seen in Fig. [121 for x=0.05 (for a 100 
GeV jet), this point is k± w 3.5GeV, for x=0.5 the point 
moves up to k± « 7GeV. 

Below this point, the behavior depends on the length, 
as seen in Fig. [Til and on the emitted gluon energy (ie 
x) as seen in Fig.fTUl For short lengths, L=lfm (with the 
gluon mean free path A=lfm as well), the higher orders 
make little difference. For longer lengths, the higher or- 
ders suppress the first order peak, and also suppress the 
small negative dip for very small k± « 0.5GeV. 

This is the effect of interference between radiation from 
multiple scattering. The first order result only has the 
interference with the creation radiation. As the length in- 
creases, the first order result pushes its collision out to an 
increasing L/2, the interference with the creation point 
becomes less severe and the result closer to incoherent. 
The interference between multiple collisions comes into 
play and keeps the result further from the incoherent. 

For a specific k±, the result increases with increasing 
L, even when weighed by 1/L as in Fig. [11] Hence, the 
length dependence of the average must be stronger than 
DC L. However, the higher orders have more of an effect 
on the result for longer lengths. The full result grows 
(slightly) more slowly than the first order result. The 
full result should approach a linear growth (ie a common 
curve when weighed by 1/ L) for very long lengths. 

The opacity series converges slowly for small k±_ and 
small X. In the appendix section [Bl the opacity series is 
seen to be oscillatory, with large successive terms almost 
canceling each other - for small x. For large x, similar 
shapes in k± are seen, just with the higher orders dra- 
matically decreasing in magnitude. In Fig. [9l even the 
fifth order result does not give a good result compared 
to the ninth order result for the shape of the small k± 
spectrum, for x=0.05. This can also be seen in Fig. [121 
where the quantum source term is plotted, removing the 
classical cascade. For x=0.05, the first order is a good 
approximation down to k± « 4GeV, summing up to fifth 
order is a good approximation down to kj_ w IGeV. For 
x=0.5, the fifth order is a good approximation all the way 
down to fcj_ ~ 0. 

For large k±, lengths L = 1 — > 4fm, and (for a lOOGeV 
jet) X > 0.05, the first order result is dominant. As k± 
is pushed down, suppression from the first order result 
is seen. We will see the effect of this is the next sec- 
tion, when plotting dN/dx, the input into an Raa{pt) 
calculation. 



2. p ^ 20GeV 

The use of a lOOGeV jet has been convenient here, 
to remove (ie push up in k±) the cuts of kinematic lim- 
its, kmax = 2x{l — x)E. x=0.05 for a lOOGeV jet is a 
gluon of 5GeV. The problematic nature of radiative en- 
ergy loss calculations for lower energy jets is illustrated 
in Fig. [I3l x=0.05 for a 20 GeV jet corresponds to a 1 
GeV emitted gluon, far from hard relative to the medium 
(remember, both the jet and the emitted gluon are con- 
sidered to be eikonal). The kmax cut in this region is so 
harsh, that we might consider the remainder of the cal- 
culation moot. While this may not be important for the 
calculation of the average energy loss, it is important for 
RAA{PT)^ which is more sensitive to the distribution in 
number of gluons, which is peaked at small x. 

A cut on a distribution is a usually a sign of pushing 
kinematic approximations too far away from their region 
of applicability. The hope is that little of the distribution 
sits in the region of kmax, so it could be varied up and 
down without changing the result. For jets at RHIC, this 
just is not possible. 



B. Results - DGLV vs GLV 

The results in the previous section were given for the 
GLV result, with the mass of the jet and the emitted 
gluon as zero. For consistency, we would hope that the 
DGLV result would give something similar for small jet 
and gluon masses. But the small mass of the emitted 
gluon has a surprisingly large effect, as shown for first 
order in Fig. ll5l and summed up to ninth order in Fig. 1161 
The result is changed well beyond an ad-hoc cut on the 
GLV spectrum at k_\_ = /i. This is worrying for the ac- 
curacy of the results. It might be necessary to include 
better considerations of the in-medium dispersion rela- 
tion of the emitted gluon. 



C. Results - dN/dx and dl /dx 

For calculations oi Raa{pt) , the important quantity is 
jet energy loss and the important distribution is dN/dx. 
For the calculation of the average energy loss, the relevant 
distribution is dl /dx — xEdN/dx. As the Gyulassy- 
Wang potential has an explicit S{q'^) factor, the energy 
of the emitted gluon equals the energy loss of the jet 
(if this factor is removed, a more considered translation 
between energy emitted and energy loss is needed, due 
to possible energy transfers in the collision). 

In order to avoid the details of a /ct > M cut, the 
results in this section use the DGLV formula, including 
a small mass for the jet and emitted gluon - despite the 
reservations expressed in the previous section. As before, 
the jet propagates through a medium of fixed density, 
fixed length, with a gluon mean free path A = Ifm and 
a Debye mass fi = O.bGeV. The cut on large k± is also 
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FIG. 1: The one-gluon number (dN/dx - left) and energy distribution (dl/dx -right) of radiation using DGLV. The upper plot 
shows the short length L=lfni for a 20 GeV jet; the lower plot shows a longer length L=4fra for a 100 GeV jet. For both, 
higher orders than those shown are negligible. 




Light quark 



L=4fm,p=20 GeV ^ 
6 



0.6 0.8 1..0 0.2 





FIG. 2: The one-gluon number (dN/dx) and energy distribution (dl/dx) for a light quark and a bottom quark jet with 
momentum 20 GeV, with L=4fm, at first order and summed up to higher orders. 
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important - here the value given by Djordjevic jT^I is 
used, k±^^ax = 2a:(l - x)E. 

The figures (Fig. [1] Fig. [2]) show that the first order in 
opacity dominates the resuhs. For lengths up to L=4fm, 
the second and third order provide small alterations on 
top of this, reducing the number of giuons emitted at 
small X. The orders above third order are negligible in 
considerations of dN/dx. 



V. CONCLUSIONS 

The opacity expansion is a Dyson expansion of a 
Schrodinger-like equation. It provides both a derivation 
of that Schrodinger equation and a method of solution 
to almost arbitrary accuracy. GLV II gave a closed form 
expression for the all orders result. This result can be 
decomposed into a quantum source term and a classical 
cascade of emitted giuons. 

A new numerical evaluation of the GLV all orders re- 
sult has been presented. Monte-Carlo integration is well 
suited to the high order integrals, and a basic version 
has been implemented to give results up to ninth order 
in opacity. These results give qualitatively the same con- 
clusions as in GLV II: for the lengths and energies of 
interest, the first order in opacity is the largest contri- 
bution. The higher orders reduce the gluon emission for 
longer lengths and small k±^ . This small k± region change 
affects giuons below k, 4 GeV. This has little effect on the 
average energy loss of jets above « 20 GeV. Raa{pt) is 
dependent on the number distribution after a multi-gluon 
emission convolution, so will be affected more than the 
average energy loss. 

Further investigation is needed to give full Raa{pt) 
results (and other observables). It is possible that the 
inclusion of higher orders may have an effect differen- 
tial in mass, as coherence is affected by the mass of the 
jet. Work is also needed on the evaluation without the 



uncorrelated medium assumption, and on the result in 
a medium of arbitrary varying density. This calculation 
is numerically more intensive, but is still possible with 
today's computing power. 

The difference between the GLV and DGLV distribu- 
tions in dN/dxdk for a light quark jet is worrying. The 
inclusion of a small gluon mass, nig « 0.3GeV, makes 
a surpisingly large difference. Further analysis of the 
medium effects on giuons with oj < 3GeV, k± < 3GeV is 
needed. 

With numerical optimization and greater computing 
time, orders in opacity beyond ninth are easily accessi- 
ble. While this is likely not important for RHIC or LHC 
applications (for perturbative input parameters such as 
A ~ Ifni), it should provide important insight into the 
phenomenon of coherent radiation in thick media, and 
possible contact with the approximations made in the 
evaluation of the BDMPS result. 

Further improvements to the formalism are necessary. 
The derivation in GLV II applies to small x, eikonal 
kinematics and no energy exchange with the medium. 
Pushing beyond these assumptions is likely necessary for 
quantitative work, especially for emitted giuons of 1-2 
GeV. The eikonal assumption applies both to the jet and 
the emitted gluon, a poor assumption in the region of 
interest. Despite its complicated derivation from QCD 
diagrams, the GLV II operator recursion from order to 
order in opacity is remarkably simple, and corresponds 
to a remarkably simple Schrodinger-like equation (as in 
BDMPS and Zakharov's work). There is a tempting 
substitution of the 2D, transverse vectors into 3 or 4D 
vectors, the longitudinal momentum shift into a general 
energy shift for an arbitrary potential. However, much 
work is needed to prove or disprove a simple result for 
the more general kinematics that are desired. 

Acknowledgments: Valuable input from Miklos Gyu- 
lassy is gratefully acknowledged, as well as discussions 
with Ivan Vitev. 
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FIG. 3: Results from the new Monte-Carlo evaluation, to be compared to fig 4 from GLV II \1^, for L=5fm, x=0.05 of a 
lOOGeV jet. The left plot shows the result for first, first plus second, and summed up to third orders; the middle and right 
show the effect of orders up to ninth order (left) and the DGLV inclusion of a jet and gluon mass (right). 
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FIG. 4: Results from the new Monte-Carlo evaluation, to be compared to fig 4 from GLV II [lO|, for L=5fm, x=0.2 of a lOOGeV 
jet. The left plot shows the result for first, first plus second, and summed up to third orders; the middle and right show the 
effect of orders up to ninth order (left) and the DGLV inclusion of a jet and gluon mass (right). 






n=1+2+3 




x=0.5 




0.6 


— ■n=1 + ...+ 6 
n=1 + ...+ 9 








0.4 










0.2 


/ 


^ \ 


V^^ 




0.0 
0.2 


-==y 




^ 




10 


100 


1000 


10 



_GLV:n=1 + ...6 
DGLV:n=1 + ...6 




FIG. 5: Results from the new Monte-Carlo evaluation, to be compared to fig 4 from GLV II [10], for L=5fm, x=0.5 of a lOOGeV 
jet. The left plot shows the result for first, first plus second, and summed up to third orders; the middle and right show the 
effect of orders up to ninth order (left) and the DGLV inclusion of a jet and gluon mass (right). 
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FIG. 6: The dN/dxdk GLV results for individual orders in opacity 1,2,3 for a 100 GeV massless quark jet. Shown is the result 
for the full result including the classical cascade, and the result only for the quantum source term. 
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FIG. 7: The dN/dxdk GLV results for individual orders in opacity 4,5,6 for a 100 GeV massless quark jet. Shown is the result 
for the full result including the classical cascade, and the result only for the quantum source term. 
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FIG. 8: The dN/dxdk GLV results for individual orders in opacity 7,8,9 for a 100 GeV massless quark jet. Shown is the result 
for the full result including the classical cascade, and the result only for the quantum source term. 
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FIG. 9: The effect of the curtailment of the opacity series. The emitted gluon is at x=0.05, ie 5 GeV. The resuhs are compared 
to the Bertsch-Gunion incoherent hmits. 
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FIG. 10: The radiation spectrum dN/dxdk± summed to different orders for different values of x. The results are weighed by 
x/L, to give converging results for large k± - the Bertsch Gunion, incoherent answer scales with L/x. 
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FIG. 11: The length dependence of the dN/dxdk± (at fixed x) radiation spectrum summed to different orders for a gluon with 
x=0.05 (ie 5 GeV). The answer is weighed by x/L, as the incoherent answer scales with L/x. The result increases with length, 
giving a length dependence stronger than oc L, but the summed result also deviates further from the first order result. 
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FIG. 12: The radiation spectrum dN/dxdk± at fixed x summed to different orders, witli (Eq. (|45p ) and witfiout (Eq. (|46p ) the 
ciassical diffusion terms. 



APPENDIX D: DN/DXDK FOR P=20 GEV 
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FIG. 13: Radiative distributions (quantum source term only, Eq. (|46|) ) for a massless quark jet of momentum 20 GeV. The 
sharp cut at small x is the kinematic cut on k±, given by Djordjevic as kmax = 2a;(l — x)E. 
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FIG. 14: The jet energy dependence of the emitted gluon distribution, for a fixed gluon energy with only the quantum source 
term, Eq. (|46|) . The distributions are very similar, except for a change in the kmax and a small effect of changing qmax = y/QTE. 
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FIG. 15: First order in opacity radiative result, for GLV {M = rUg = 0) and DGLV(Mjet = A'/2 = 0.25 GeV, nig = /i/V2=0.354 
GeV), for the quantum source term only , Eq. (|46p . The inclusion of the masses has a surprisingly large effect for gluon energies 
and fexs up to ~ 3GeV. 
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FIG. 16: Similar to the previous plot (GLV vs DGLV comparison Fig. I15p . but for the result summed up to ninth order in 
opacity (again, quantum source term only, Eq. (|46|) ). 
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